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We introduce a multivariate diffusion model that is able to price derivative securities 
featuring multiple underlying assets. Each underlying shows a volatility smile and is 
modeled according to a density-mixture dynamical model while the same property holds 
for the multivariate process of all assets, whose density is a mixture of multivariate ba- 
p I \ sic densities. This allows to reconcile single name and index/basket volatility smiles in 

a consistent framework. Rather than simply correlating one-dimensional local volatility 

•i-H \ 

M-H models for each asset, our approach could be dubbed a multidimensional local volatility 

approach with state dependent diffusion matrix. The model is quite tractable, leading to 
a complete market and not requiring Fourier techniques, contrary to multivariate stochas- 
tic volatility models such as Wishart. We provide a semi-analytic formula for the price 
of European options on a basket /index of securities. A comparison with the standard 
' approach consisting in using Monte Carlo simulation that samples simply-correlated suit- 

ably discretized one-dimensional paths is made. Our results show that our approach is 
promising in terms of basket option pricing. We also introduce a multivariate uncer- 
£f} • tain volatility model of which our multivariate local volatilities model is a multivariate 

markovian projection and analyze the dependence structure induced by our multivariate 
dynamics in detail. A few numerical examples on simple contracts conclude the paper. 
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1 Introduction 

It has been known for a long time that the Black-Scholes geometric Brownian motion 
model does not price all European options quoted on a given market in a consistent 
way. In fact, this model lies on the fundamental assumption that the asset price volatility is a 
constant. In reality, the implied volatility, namely the volatility parameter that, when plugged 
into the Black-Scholes formula, allows to reproduce the market price of an option, generally 
shows a dependence on both the option maturity and strike. If there were no dependence on 
strike one could extend the model in a straightforward fashion by allowing a deterministic 
dependence of the underlying's instantaneous volatility on time, so that the dynamics could 
be represented by the following stochastic differential equation (SDE): 



at being the deterministic instantaneous volatility referred to above. In that case, recon- 
struction of the time dependence of at would follow by considering that, if v(Ti) denotes the 
implied volatility for options maturing at time Tj, then 



Implied volatility however does indeed show a strike dependence; in the common jargon, 
this behavior is described with the term smile whenever volatility has a minimum around 
the forward asset price level, or skew when low-strike implied volatilities are higher than 
high-strike ones. In the following we will loosely speak of both effects as "volatility smile". 

In recent years, many researches have tried to incorporate the smile effect into a consistent 
theory. Generally speaking four main streams of investigation can be identified in a univariate 
setting. 

A first approach is based on assuming an alternative explicit dynamics for the asset-price 
process that by construction ensures the existence of volatility smiles or skews. An example 
is the CEV process proposed by Cox |17| and Cox and Ross [18j. A different example is 
the displace diffusion model by Rubinstein |45j . A general class of problems was presented 
by Carr, Tari and Zariphopoulou [16j. In general the alternative explicit dynamics does not 
reproduce accurately enough the market volatility structures, since it is based on quite stylized 
dynamics. 

A second approach is based on the assumption of a continuum of traded strikes [1] . This 
was extended yielding an explicit expression for the Black-Scholes implied volatility as a 
function of strike and maturity |20} [2T| [22} 123] , This approach however needs a smooth 
interpolation of option prices between consecutive traded strikes and maturities. Explicit 
expressions for the risk-neutral stock price dynamics were also derived by minimizing the 
relative entropy to a prior distribution [lj and by assuming an analytical function describing 
the volatility surface [14]. 



dS t = ixStdt + a t S t dWt 
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Another approach is an incomplete market approach, and includes stochastic volatility 
models |33[ [34"1 |4"6] and jump-diffusion models |44j . 

A further approach consists of finding the risk-neutral distribution on a lattice model for 
the underlying that leads to a best fit of the market option prices subject to a smoothness 
criterion |13[ I36| . This approach has the drawback of being fully numerical. 

A number of the above approaches is described for the foreign exchange market in Lipton 
|40| . see also Gatheral }27] who deals further with volatility surfaces parameterization. Recent 
literature also focused on both short time and long time asymptotics for volatility models: 
we just cite |28| as a reference for small time asymptotics for local volatility models, and |26| 
for large maturities asymptotics in the well known Heston stochastic volatility model, while 
pointing out that the volatility asymptotics literature is much broader. 

In general the problem of finding a risk-neutral distribution that consistently prices all 
quoted options is largely undetermined. A possible solution is given by assuming a particular 
parametric risk-neutral distribution dependent on several, possibly time-dependent, parame- 
ters and use the latter in conjunction with a calibration procedure to the market option prices. 
In a number of papers, Brigo, Mercurio, Rapisarda and Sartorelli [SJ [U [71 El [10] proposed a 
family of models that carry on dynamics leading to a parametric risk-neutral distribution 
flexible enough for practical purposes. It is relatively straightforward to postulate a mixture 
distribution at a given point in time, but it is less so finding a stochastic process that is con- 
sistent with such distribution and whose stochastic differential equation has a unique strong 
solution. This is the approach adopted by the above papers. This family of models is sum- 
marized for example in Musiela and Rutkowski |42| . or Fengler [25], see also Gatheral |27j . 
Formally, this is part of the alternative explicit dynamics branch of models but is typically 
much richer than the models listed above, leading to a practically exact fit of the volatility 
smile while retaining analytical tractability. 

The aim of this paper is to incorporate the effect of the volatility smile observed on the 
market when pricing and hedging multiasset securities, while retaining sensible single asset 
volatility structures. A whole lot of such structured securities is nowadays offered to insti- 
tutional and retail investors, in the form of options on baskets of stocks/FX rates and on 
combinations of forward interest rates such as e.g. European/Bermudan swaptions. In our 
approach we remain within a local volatility model for the individual assets composing the 
underlying of the option (be it a basket of stocks or a swap rate) that has proved to be quite 
effective in accounting for the observed single assets' smiles, but we move one step beyond 
the approach used on the street when writing the joint multi- asset dynamics. Given uni- 
variate local volatility (one dimensional diffusion-) models for each asset, a basic approach is 
introducing instantaneous correlations across the brownian shocks of each asset, leading to 
what we call the Simply Correlated Mixture Dynamics (SCMD). For practical implementa- 
tion, one would then discretise the one-dimensional single asset SDEs through, say, Euler or 
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higher order numerical schemes |38j . feeding correlated instantaneous Browinan shocks into 
the scheme. In this paper we adopt a different approach and we incorporate the correlation 
in a new scheme that enjoys analytic multivariate densities and a fully analytic multivariate 
dynamics through a state dependent non-diagonal diffusion matrix. In so doing we are able 
to sample a new manifold of instantaneous covariance structures (and a new manifold of dy- 
namics) which ensures full compatibility with the individual volatility smiles and overcomes 
the difficult problems created by the lack of closed form formulas for prices and sensitivities. 
We call the resulting model Multi Variate Mixture Dynamics (MVMD). 

The level of tractability in MVMD for both single assets and indices/baskets is much 
higher than with multivariate stochastic volatility models such as Wishart models, for which 
we refer for example to |29[ I19j and references therein. Furthermore, the MVMD model leads 
to a complete market and hedging is much simpler. It is practically a tractable and flexible 
multivariate local volatility model that has the potential to consistently calibrate univariate 
and index volatility smiles through a rich but at the same time transparent parameterization 
of the dynamics. 

The traditional approach for pricing European style derivatives on a basket of the multi- 
dimensional underlying, in a SCMD type model, uses a Monte Carlo method that can be very 
slow as it involves intensive time discretization, given that correlation can only be introduce 
at local shocks level. With this paper we fill this substantial gap in option pricing and pro- 
vide, with MVMD type models, a semi-analytic solution to the option pricing problem where 
the price can be quickly and accurately evaluated, something that practitioners value greatly, 
especially in Risk Management area. 

In order to develop a feeling for the performance of our approach, we test it on a few cases. 
We distinguish two types of basket options: An option on a weighted arithmetic average of a 
basket and an option on a weighted geometric average of a basket. We also consider spread 
options with a smile structure of implied volatilities. We compare the prices generated by 
MVMD to those obtained by the SCMD model with analogous parameters. 

We then introduce a Multivariate Uncertain Volatility Model (MUVM). We show that the 
MVMD model is a markovian projection of the MUVM. MUVM thus gives the same European 
option prices as MVMD and can be used instead of MVMD to price European options also 
in the multivariate setting . MUVM features the same dependence structure as the MVMD 
model. The related copula is a mixture of multivariate copulas that are each a standardized 
multivariate normal distribution with an appropriate correlation matrix and margins. 

The paper is organised as follows. In Section [21 we present a brief review of the approach 
to single asset smile modeling that has been developed in [5j [TJ [8] . In Section [3j we provide 
examples of typical securities that need a multivariate setting for proper pricing. Section 
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[4] considers the extension of the single asset model to the multivariate framework with a 
thorough discussion of the implications for the dynamics stemming from a naive approach 
(SCMD) and from ours (MVMD). This section also provides a semi-analytic formula for 
pricing European options on a basket and spread options under MVMD. In Section we 
introduce a new model that we call "Multivariate Uncertain Volatility Model" so that our 
model is a multivariate markovian projection of it. In Section [6l we derive the copula function 
associated to our model. In Section [71 we illustrate the results of pricing European option on a 
weighted arithmetic average of the underlyings with positive weights, European spread option 
and European option on weighted geometric average in MVMD and SCMD frameworks and 
we compare the results. Conclusions and suggestions for future research are given in the final 
section. 

2 The Mixture Dynamics (MD) Model 

We fix a time T > and denote by P(0, T) the price at time of the zero-coupon 
bond maturing at T. Let (Q,J-, P) be a probability space with a filtration (J~t)te[o,T] that 
is P-complete and satisfying the usual conditions. We assume the existence of a measure 
Q equivalent to P called the risk neutral or pricing measure, ensuring arbitrage freedom in 
the classical setup, for example, of Harrison, Kreps and Pliska |31[ I32| . At times, it will be 
convenient to use the T-forward risk-adjusted measure Q T rather than Q. 

The MD model is based on the hypothesis that the dynamics of the asset underlying a 
given option market takes the form 

dS(t) = n(t)S(t)dt + u(t, S(t))S(t)dW(t) (2.1) 

under Q with initial value So- Here, \x is a deterministic time function, W is a standard Q 
Brownian motion and v (the "local volatility") is a well behaved deterministic function. In 
order to guarantee the existence of a unique strong solution to the above SDE, v is assumed 
to satisfy the linear growth condition 

u 2 (t, x)x 2 < L{1 + x 2 ) uniformly in t (2.2) 

for a suitable positive constant L. 

Consider TV purely instrumental diffusion processes Y l (t) with dynamics 

dY\t) = ^(t)Y l (t)dt + Vi(t, Y^t^Y^tjdWit) (2.3) 

with initial value Y l (0), marginal densities p\ and V{ satisfying linear growth conditions, where 
each r*(0) is set to 5(0). 

Remark 1 The reader should not interpret the Y % as real assets. They are just instrumental 
processes that will be used to define mixtures of densities with desirable properties. 
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The marginal density pt of S(t) is assumed to be representable as the superposition of the 
instrumental processes densities p\ [6j [TJ [8] : 



Pt 



\p\ with Aj > 0,Vz and \ = 1. 



(2.4) 



The problem of characterizing v can then be cast in the following form: is there a local 
volatility v for Eq. (|2,ip such that Eq. (|2,4p holds? Purely formal manipulation of the related 
Kolmogorov forward equation 



dpt d 1 d 2 2 >2 \ 

~at + ^ xpt) ~ 2M {u {t > x)x Pt) = 

and of analogous equations for the p\'s shows that a candidate v is 



Eli^vf(t,x)pl(x) 



v(t, x) 



(2.5) 



(2.6) 



A digression is in order now: suppose that everything in the above approach works fine, 
and with a judicious choice of the Vi's the unique strong solution for the dynamics of S exists. 
Then the pricing of European options on S is immediate: Let O be the value at t = of 
an European option with strike K and maturity T. O is given by O = ^i=i ^iOf, where Oj 
is the European price associated to a hypothetical instrumental dynamics (|2.3p . The option 
price O can be viewed as the weighted average of the European option prices written on the 
processes Y % . Due to the linearity of the derivative operator, the same convex combination 
applies to all option Greeks. As a consequence, if the basic densities p\ are chosen so that the 
prices Oi are computed analytically, one finds an analytically tractable model. 

The most natural choice for the (Y l ,Vi,pl) triplet is : 
' y*(0) = 5(0) 



Vi(t,x) = o-j(i) 



V; 



(2.7) 



p\{x) 



2irxV t (t) 



exp 



2Vi(t) 



In 



with (jj deterministic (lognormal mixture dynamics, LMD). 

Brigo and Mercurio [7] proved that, with the above choice and additional nonstringent 
assumptions on the Oj, the corresponding dynamics for St admits indeed a unique strong 
solution. A greater flexibility can also be achieved by shifting the auxiliary processes' density 
by a carefully chosen deterministic function of time (still preserving risk-neutrality). This is 
the so-called shifted lognormal mixture dynamics model [8]. 
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The above description gives a sufficient basis for presenting our generalization of the LMD 
to the multivariate setting, as before at first on the basis of pure formal manipulations, and 
then with full rigor, with the specific aim of finding a method to infer the "implied volatility" 
of a basket of securities from the individual components and/or an explicit dynamics for the 
multi-asset system. Later in the paper, formal proofs of the general consistency of the model 
and of the existence and uniqueness of the solution to the multivariate version of Eq. (|2.ip 
will be provided. 

3 Options on Baskets: Motivating multivariate models 

A generalization of LMD to the multivariate setting aims to be able to compute the smile 
effect on the implied volatilities for exotic options depending on more than one asset, such as 
a basket options. Clearly, analogous techniques apply to indices. 

3.1 Basket option 

A basket option is an option whose payoff is linked to a portfolio or "basket" of underlying 
assets. We can distinguish two types of basket option: 

• An option of weighted arithmetic average of the under lyings: 

B t = J2^kS k (t), (3.1) 
k 

• An option of weighted geometric average of the under lyings: 

n 

Br- l[S k (t) w * 

,fc=i 



(3.2) 



where Sk is the k-th component of the basket. Typically the basket is consisting of several 
stocks, indices or currencies. Less frequently, interest rates are also possible (Sk could rep- 
resent a forward rate process Fk in the Libor Market Model (LMM) and instead of (|3.ip we 
could have a more complicated expression representing a swap rate). 

Such options have the most varied nature: from the plain European call/put options on 
the value of the basket at maturity T, to options somewhat more complicated, such as Asian 
options on the basket, Himalaya options, rainbow options and so on. 

The weights (u>k)k i n (|3.ip can be negative. When the basket (|3.ip contains short positions 
it is called spread and the option known as a spread option is written on the difference of 
underling assets. The weights (wi~)k=i ...n i n (|3 . 2|) are positive. 

It is instructive to view a basket option as a standard derivative on the underlying instru- 
ment whose value at time t is the basket Bt so defined. In this paper we will be concerned 
only with European options. 
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3.2 European options pricing 

Let us assume that interest rates are constant and equal to r > . We also assume the 
existence and uniqueness of a risk-neutral pricing measure Q that is equivalent to P under 
which discounted asset prices are martingales, implying the absence of arbitrage (Q is also 
equal to Q T as interest rates are constants). 

According to the Black-Scholes pricing paradigm |31| [32] . the price II of an European 
option at initial time t = is given by the risk-neutral expectation: 



where the exponential factor takes care of the discounting and uj = ±1 for a call/put 
respectively. Bt is the underlying instrument (can represent the value of the basket) at 
maturity T, K is the strike. 

The fundamental difficulty in pricing basket options on a weighted arithmetic average of a 
basket is to determine the distribution of the sum of underlying asset prices. Let us consider 
the basket of securities of Eq. (|3.ip . Several approximation methods have been proposed for 
options on it. Usually the basket value (|3.ip is approximated by the lognormal distribution. 
Recall that here we consider baskets with possibly negative weights, such as spreads. Hence, 
we cannot approximate the distribution of Bt by a lognormal distribution, since such a basket 
can have negative values or negative skewness. However, Borovkova, Permana and weide |12| 
show that a more general three-parameter family of lognormal distributions: shifted, negative 
and negative shifted lognormal, can be used to approximate the distribution of a general bas- 
ket. The shifted lognormal distribution is obtained by shifting the regular lognormal density 
by a fixed amount along the x-axis, and the negative lognormal - by reflecting the lognormal 
density across the y-axis. The negative shifted lognormal distribution is the combination of 
the negative and the shifted one. Note that this family of distributions is flexible enough 
to incorporate negative values and negative skewness: something that the regular lognormal 
distribution is unable to do. However, by using these approximations we do not take into 
account the internal composition of the basket value in terms of underlying assets having 
each its own dynamics. This approach structurally cannot take into account any smile effect 
on the individual under lyings' volatility, and therefore on the "basket volatility". 

In the following we will tackle the problem in a rigorous way, through the generalization 
of a dynamical model that has proven to perform quite well on some markets [H [B] and 
that is under extension to the equity markets case. 




(3.3) 
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4 An extension of the MD model to handle multivariate prob- 
lems 

To fix ideas, suppose we are faced with the following problem: we want to price an option 
maturing at T on the basket of securities given by Eq. (|3,ip or Eq. (|3,2p . Each of these 
securities will have a "smiley" volatility structure, and we strongly suspect that the basket 
will show a smile in its implied volatility, too. 

Through Eqs. (|2.1H2.7p we now have a piece of machinery that allows us to calibrate an MD 
to the implied volatility smile structure of the individual component of the basket. Suppose 
we have already calibrated the individual MDs to such smile surfaces, thus finding the local 
volatilities governing the dynamics of each Sk ■ We denote by , . . . , Yir the instrumental 
processes for asset Namely, for each asset in the basket we have a family of instrumental 
processes like (|2.3p that refer to that specific asset mixture distribution. Most often, a good 
way to price a plain European option depending on more than one asset is to use a Monte 
Carlo simulation that samples suitably discretized paths according to the drift rate of each 
component (risk-free minus dividend yield) and to the diffusion matrix given by the local 
volatility function deriving from the mixture of densities model. Therefore, assuming to have 
an exogenously computed structure of instantaneous correlations pij (computed e.g. through 
historical analysis or implied by market instruments and supposed constant over time) among 
the assets' returns, we could apply a naive Euler Monte Carlo scheme and simulate the joint 
evolution of the assets through a suitably discretized time grid T\ = • • • tjv = T with a 
covariance matrix whose component over the (t/,t/ + i) propagation interval (constant in 
the Euler scheme) is 



Cij(S*j, Sj, t) 



Efe=l X i^H T l)PY^( n )( S i) 



Ek=i X j a j( T i)PYH n )( s i, 



Pir (4-1) 



It can easily be shown that the approach of Eq. (|4,ip called SCMD is consistent with both 
the individual dynamics induced by a MD model for each underlying asset and with the 
imposed instantaneous correlation structure p^. However, besides the practical possibility of 
controlling the instantaneous correlation, and that the number of base univariate densities to 
mix does not increase with the number of underlying assets, one must be aware of its main 
limitations: 

• By imposing that the instantaneous covariance of the multidimensional process be of the 
form lj[4.1\ ) one is moving only within a given manifold of the possible local volatility 
structures for the multidimensional process. 



• The fact that we are considering options of the European type implies that we do not 
really need the full dynamics when it comes to actually computing the price of a basket 
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option. Indeed, the payout at maturity depends only upon the values of the assets at 
time T, i.e., upon the values Sk(T), \/k, regardless of the history of prices. So in order 
to compute the risk-neutral expectation in h3. 3\) giving the price II, the only information 
we need is the joint density of the process (S\(T), S%(T), . . .) of random variables under 
that particular risk-neutral measure. This density is usually called the state price density. 
But here we do not know it, so we have to generate samples from the entire path Bt for 
< t < T . The discretization time steps Ti + \ — t\ should be chosen carefully to be sure 
that the numerical scheme used to generate the discrete samples produces reasonable 
approximations. Notice that when the maturity T increases, more time steps are needed. 

One could try to do something different and approach the problem so that, under suitable 
assumptions, the individual MD models (one for each underlying asset, separately calibrated 
each on its volatility surface) could be merged so as to provide a coherent multiasset model 
that allows for a degree of (semi) analytic tractability comparable to the one typical of the 
univariate case. 

for this purpose we need to briefly revise some properties of the multidimensional ltd 
processes, and extend the existing MD Model. 

4.1 The multidimensional Kolmogorov equation 

Consider an n-dimensional stochastic process S_{t) = [S\ (t), ■ ■ ■ ,S n (t)] T (where T repre- 
sents the transpose operator) whose generic i— th component follows the SDE 

dSi(t) = iXiSi{t)dt + Si(t)Ci(t,S)dW(t) (4.2) 

where fa is a constant, W_ = \W\,--- 1 W n \ r is a standard n-dimensional Brownian motion 
and Ci(t, S) is a row vector whose components are deterministic functions of time and of the 
state of the process S_. 

Denote C{j(t,S_) = Ci(t,S_) Cj (t,S_). The associated Kolmogorov forward PDE to be 
satisfied by the probability density Pstt) OI the stochastic process S_ is 

9ps(t) , A d 1 <A d 2 r ^ 

-^T + Z^Q^^SSd ~ 5 L WS( t) ] = (4.3) 

where all functions are evaluated at (t,x) Vi > 0,x £ M n . 
With these notations S_ is given by the SDE 

dS(t) = diag{p)S(t)dt + diag(S(t))C(t,S(t))dW{t) (4.4) 

where C is an n x n matrix whose % th row is Cj. 

C must be chosen so as to grant a unique strong solution to the SDE (|4.4p . In particular, 
C is assumed to satisfy, for a suitable positive constant K, the generalized linear growth 
conditions 

l|C(i,x)|| 2 ||x|| 2 < K(l + \\x\\ 2 ). (4.5) 
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The symbol |||| denotes here vector and matrix norms. 



Consider an n dimensional stochastic process XM whose generic i th component follows 
the dynamics 

dx\ k) {t) = HiX\ k } (t)dt + x\ k) {t)af ] {t, X^)dW(t) (4-6) 

with 5?° an 1 x n matrix satisfying particular conditions ensuring that the resulting 
SDE giving the dynamic of JC. has a unique strong solution. 



Denote a {k \t, 1^) = af } (t,X (&) ) of \t, X^) T and pf ] the probability density function 



Mi 



of X^ . The associated Kolmogorov equation to be satisfied by pf' is 

[a>y(t,x)xiXjpl k \x)] = 0. 



dpt^ (x) 



Of 



d 

+ ^-^-[^ipf\x 



8=1 



2 ^— ' dxidxj 

i,3=l J 



(4.7) 



Inspired by the univariate approach which gave rise to the MD model, let us postulate 
that the density at any time t of the multivariate process S be equal to a weighted average of 



the p\ 



(k) 



N 



N 



k=l k=l 

The condition that Ps(t) satisfy Eq. ([4.3p and that each p^ satisfy the equation (|4.7p leads 
through standard algebra to the PDE 

' N N \ 



1 



& 1 



2 ^— ' dxidxj 

1,3=1 



Cij(t,X_) ^ X kPt (x) - ^ (t,x)p\ ''(x) XiXj 



k=l 



k=l 



0. 



(4.9) 



This is a rather complex PDE, that we solve by setting, for all i 



Mi 



M 



(*), 



Cij(t, x) 



T,k=l X kPt (%) 

4.2 The connection between univariate and multivariate MD 



(4.10) 



We make a strong assumption, namely that the volatility coefficient for the fe-th "base" 

(k) 

density p\ of Eq. (14 .7p is a deterministic function of time, independent of the state, and of 
the particular form a^ k \t, x) = o^(t) af\t) T . 

Under this hypothesis we already know the dynamics corresponding to Eq. (|4.7p . and we 

(k) 

can explicitly write the density p t statisfying it at all times: 



p {k \x) 



1 



(2 7 r)? v / detHW(t)n™ =1 x i 



exp 



Z T {^ k \t))- l x 



(4.11) 



1 \t can be easily proven that the most general solution of the equation £V dx d gx fij (x) = has a Fourier 
transform satisfying (q, /(q)q) = 0. The only matrix function /(q) satisfying it and infinitely differentiable 
with respect to q is constant. This constant must be zero in order to have finite first and second moments of 
the multivariate density Ps{t)- 
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where 3^ k \t) is the n x n integrated covariance matrix of returns for the many components 
of the process 

Sg°(*)= f<??\ S )<?f\s) T ds (4.12) 



o 



^w(fc) i s assumed to be invertible at all times and correlation is included into the vector 
components) and 

f* a [k)2 (s) 

Xi = lnxi-]nxi(0)-Hit+ 1 v ' ds. (4.13) 



Life is not simple unless one is ready to make the further assumption o : \ k \t)W^ k \t) T = 
\\®i k \t)\\ \\^ k \^)\\pij- The f a ct that the densities will get mixed up through Eq. (|4.8|) will 
have dramatic consequences on the actual structure of correlations, both instantaneous and 
average. But first, let us prove that under a further assumption we can be fully consistent 
with the dynamics specified by the MD model for the individual assets. 

Let's assume that we have calibrated an MD model for each S{(t): if PsJt) is the- density 
of Si, we write 

psivto =E A ^ fe (*)( a; )' with A * ^ °> yk and E A * = 1 ( 4 - 14 ) 

k=l k 

where K 1 , K ' are instrumental processes for Si evolving lognormally according to the 
stochastic differential equation: 

dYf(t) = mYf(t)dt + af(t)Yf(t)dZi(t), (dZi, dZj) = Pijdt (4.15) 

with density Pyk^y 

For notational simplicity we will assume that the number of base densities Ni will be the 
same, N, for all assets. The exogenous correlation structure pij is given by the symmetric, 
positive-definite matrix R. 

The most natural tentative choice for the base densities of Eq. (|4.8j) is 

N 

*>«„(*)= E »f-*Vw ^«')1 T<£) (4 ' 16) 

k\,k2,---k n =l L J 



with 



Pi ki , Mx) = (4.17) 



1 

exp 



;(Ja-*B)Tn(fa-*ii}m-lj(*i-*ri) 



(2tt)? v /dets(fcl '" fcn) W n ?=i^ 
Here, H^ 1 '"^™)^) is the integrated covariance matrix whose element is 

E ^) {t)= j\k i{s)a k j{s)pijds (41g) 
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and, generalizing Eq. (|4.13p 

~(ki-k n) =lnx ._ lnx .( 0) _ w t + / ^f^ds. (4.19) 



~(kl-*») • ■ / K («) 



^ 

Then, from (|4,10p . the multivariate extension of MD model that we call Multi Variate 
Mixture Dynamics (MVMD) model is given by Eq. (|4.4p where C is defined by 

Et...,*„=i AN-Afr V^.-^(t) P r fe fe) 
C(t,,)C^,*) = -g— j- 1 1 " J (4.20) 

Vi, i, _i A/. ..A n p r fe , , iT X 

^k u ...,k n -i i I (*).-•• >>^f ™ C*) J 

where 

y*i *» (t)= [^(t)^^(t)].. (4.21) 

Putting notational complexity aside, what we ultimately did is to mix in all possible ways 
the component densities for the individual assets, still ensuring consistency with the starting 
models for the components assets, and imposing the instantaneous correlation structure R at 
the level of the constituent densities. 

This is indeed consistent with the mixture of densities models through which we have 
specified the dynamics of the single components of S_: 

Proposition 2 For any test function f : R — > R and any t > 0, under the assumptions of 
Eqs. l[4-16\ )- (4-19 ), the expectation of f(Si(t)) is 



Mf(Si(t))} = I dx t f( Xi ) Psdt) (x t ) (4.22) 

PSi(t) bein 9 9^en by Eq. 



Proof. The proof is trivial: it is enough to compute the multiple integral 
^o{f(Si{t))} = J dxi- • • / dxi- • • / dx r J(xi) P s(t)(x) 

(4.23) 

= Ek u k2,-k n =l A l' ■ ■ ■ fdx 1 ---fdx i ---f dx n f{x i )p^ l ^ ^ ,t(x) 

Integrating out all variables but x% in each of the integrals in the right hand side we have 



(4.24) 



E^=i >(* J dxif(xi)p Y k i{t) (xi) =fdxif(xi)p m (xi) 



since by the condition that probability integrate up to one we know that Y2k=l = 1 for all 
i. m 

An immediate consequence of this is the following 
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Corollary 3 Under the assumptions of Eqs. l[4-16\ )- ft4-19ty , all moments, at any order, at 
any time, of the individual components are exactly reproduced. 

In other words the theory we have built is fully consistent from the statistical viewpoint 
with the MD model postulated at the beginning for each component at any time. This implies 
also dynamical consistency, which is also directly proven in the following 

Proposition 4 Under the assumptions of Eqs. fa4-16\ )- (4.1ty , the unconditional dynamics of 
each component of S_ is exactly the same as in the corresponding MD model. 

Proof. Take for instance the first component: integrate Eq. (|4.7|) with respect to x 2 , ■ ■ ■ x n : 



Wt [fdx 2 ---fdx nP r H 



Y^(t),...,Yn n (t) 



T + 



dxi 



+ ••• + / dx n 



„ 2 x 2 fdx 3 ...fdx n p^ kl{t) ^ ^ 



1 d 2 



dx n 



o-* 1 (t) 2 xlJdx 2 ---Jdx nP r ki 



»nXnfdx 2 ...fdx n „ lP ^ kl{t) ^ (t)] : 



Y^{t),...,Y^{t)] 



9 2 



18x2 

j'r/.r, J\/.r,— • 



l 1 (*)o"2 2 (t)pi 2 xix 2 fdx 3 ---f dx n p r . 



Y^(t},...,Y^(t)]' 



°2 2 (t)^ 3 (t)P23X 2 X 3 f dX4 ■ ■ ■ / dx n p r . 



Note ^\ Xl ) = fdx 2 ...dx nP r ki{t) ^ Y>{t) 



dx 



2 • • • 



Y^(t),...,Yn n (t) 

t(x) and denote 
dXi-i I dxi+i ■ ■ ■ dx n Pr k , i 



+ 



T (X) 



-- 

(4.25) 
(4.26) 



the above equation becomes 

d A lk i) , d 



+/"2 



^ dxi 



X 2 TT- 



(lfcl,2fc 2 ) 



X2=00 



x 2 =0 



H VHr. 



(lfcl,nfe n ) 



1 d 2 
2dx\ 

a 

dxi 

d 
dxi 



x„=0 



a'l 1 (t)o-2 2 (i)pi2^ia:2vrf fcl ' 2A:2) 



(4.27) 



(7f 1 (t)a*»(t)^ n xi a ; n 7r t (1 * 1 ' n *" ) 



x z=° 

x n =co 



o- 2 2 (t)o- 3 3 (t)p 23 x 2 x 3 n ( t 



(2fc 2 ,3fc 3 ) 



£n=0 
X2=X 3 =00 

X 2 =X 3 =0 



+ 







All terms evaluated in zero and infinity in the above equation vanish because in order to 
have finite first and second moment with respect to all components, the following limits hold 



linx^+oo XiPs(t) (x) = 
lim^+oc xfps( t ){x) = 



lim^o XiPs(t) {x) = 
lirn^o xfpsit) (x) = 



(4.28) 
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and similar conditions hold for the derivatives of Ps(t)- 

Each term in Eq. (|4.7j) for each base multivariate density therefore reduces to the corre- 
sponding term of the Kolmogorov equation for the one-body density of x\. ■ 

So far we have only provided a candidate density and dynamics for the multivariate prob- 
lem. We must still prove that our choice does indeed lead to a unique strong solution of the 
vector SDE for S. Moreover a deep analysis of the consequences of this choice has to be 
performed. 

Existence and uniqueness of a strong solution follows from the usual technical assumptions 
for the individual volatilities af': supposing they are all bounded from above and below [6], 
standard algebra yields 

' < IICII 2 = Et j= iC^t,S_r < n*d\ 

< a = 'mf t > {mm i=1 ... nM=1: ... N }a^(t)} (4.29) 

a = sup 4 > {max i= i... niA ; != i...Ar}crf l (t)} 
and the linear growth condition holds. 

4.3 Dimensionality issues 

The computational scheme shown above ensures full consistency between the single-asset 
and the multi-asset formulations of the mixture of lognormal densities' model. It must be 
borne in mind, however, that the number of "base" multivariate densities of the formulation 
of Eqs. (|4.16p - (|4.19p explodes as N n if we have N base univariate densities for each of the 
n underlying assets (more generally, if asset i relies on a single-asset mixture theory based 
on Ni densities, the number of multivariate densities entering the superposition amounts to 
nr=i^)- This "combinatorial explosion" seems to limit the applicability of the theory to 
baskets made of very few assets. 

However, as already observed elsewhere [8] two empirical facts appear in the univariate 
mixture of densities model, that encourage the application of the model to real world multi- 
variate settings. They are briefly summed up here: 

• the number of base densities N needed to reproduce accurately enough the implied 
volatility surface for a single asset is typically 2 to 3; 

• there appears to be a clear hierarchy between densities composing the mixture, dictated 
by the weights borne by each density in the superposition: in fact, typically one 
density takes up most of the weight, the second takes up most of the remaining weight 
(remember that X^fcLi ^fe = 1) an d t ne l as t weighs little compared to the first two. 
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The consequences of the first issue are evident: the base in the power law N n is of the 
order of two/three. This is not enough to completely solve the explosion problem: taking 
N = 3 and n = 8 still implies that in order to compute the price of an European option on 
the basket, we should compute 6561 multidimensional integrals. 

However, the second point ensures that most of the multivariate coefficients A^ 1 • • • A^ n 
result from the product of the smallest A, thus rendering the corresponding terms in the 
expansion of Eq. (|4.16p negligible. Given any < K < 1, a possible solution can therefore be 
to approximate Eq. (|4.16p through 

N ( n \ 

fe 1 ,fc 2 ,-fen=l \i=l / L 1 " J 

such that n™=i A j J>K 

K therefore plays the role of a "cutoff parameter" that ensures that only significant contri- 
butions to the multivariate expansion are retained; in order to preserve normalization of the 
resulting density, 

n yin ykj 

11 ? = -■ (4.3D 

E II >■/ 

k\,k2,---k n =l 

such that n?=i -V >K 

Note that Ps(t)(%) = Ps(t){%i K = 0)> whereas increasing k decreases the number of base 
multivariate densities in the approximate expansion of Eq. ([4.3UP : K therefore controls the 
tradeoff between the accuracy in the approximation and the computational efficiency. 

In order to have an estimate of the computational gain due to a choice of K ^ 0, we can 
compute how the volume in n-dimensional space of the region k < ]Xi=i x i ^ 1 scales for 
fixed cutoff as a function of n > 1: the recursive law is 

V n (K)= f dx x f dx 2 --- f dx n = Vn-^K) + K(lnK)"- 1 (4.32) 

Jk J^- J n-1! 

setting Vq = 1 conventionally. 

Now, let us neglect the striking feature that there exists a strong hierarchy between com- 
ponents in the univariate mixtures of densities (point two above). Suppose instead that we 
are in a less favorable case, namely that the density of coefficients A of the mixture model for 
each asset is uniform and equal to p {i.e. the distance on the [0, 1] interval between consec- 
utive A is equal to =); then, the density of coefficients in the multivariate theory is p n . An 
estimate of the number of multivariate densities involved in the expansion of Eq. (|4.30p is 
9tn(ft) = Vn{ K )p n ■ To give an example, if k = 5% and p = 3, the number of densities has a 
maximum at 9T n (5%) ~ 80 for n ~ 8: neglecting the densities that contribute to 5% of the 
normalization already yields much less than the full 6581 set of eight variate densities. 
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4.4 European options pricing 

Suppose that S_ represents the vector of underlying asset prices composing the underly- 
ing B in Eq. (|3.ip or Eq. f|3 . 2[) . A conventional scheme for pricing a plain option on the 
underlying basket in a way fully consistent with individual local volatilities would require, 
according to a SCMD type approach, a sufficiently fine time discretization coupled with a 
Monte Carlo integration with instantaneous covariance given by Eq. (|4.ip (or by more com- 
plicated discretization schemes for SDEs, see e.g. Milstein's |38j). Our MVMD model allows 
instead to compute the option price (|3,3p through a set of single-step Monte Carlo integrations 
(one integration for each combination (fcj., ■ ■ ■ ,k n )). Indeed since the terminal distribution of 
S_(T) is known, the MVMD model allows to evaluate simple claims on a basket without time 
discretization. Thus, using the MVMD approach, one can reduce the computational time 
significantly. But the actual consequences of this approach are wider, in that they affect the 
many-body dynamics in a deeper way. 

Remembering (|4.16p . it is straightforward to obtain the model option prices in terms of the 
option prices associated to the instrumental processes (momentarily thought of as underlying 
assets) (l^ fc )i = i > ... >n ,k=l,- ,N- 

4.4.1 Option on a weighted arithmetic average of a basket 

Let us consider an option of European type on the basket of securities of Eq. (|3.1|) with 
maturity T and strike K. The discount rate is denoted r. It is supposed to be constant. Then, 
if oj = 1 for a call and ui = — 1 for a put, the option value (|3.3p can be written as 

+ 

Ps(T)(xx,...,x n )dx 1 ...dx n (4.33) 

where p$(T) is the joint density of the random variables Si(T),. . ., S n (T) and is given by Eqs. 
(T4~16|) - (j4~T9]) . Finally we have 

N 

n= X l 1 --- X n *l,..,fcn (4.34) 

k\ , . . . , kn — 1 

where &k 1: ...,k n denotes the European option price associated to the basket Y17=i w i^i l - 

When the value of the basket (|3.ip contains short positions we are dealing with spread 
options. 

4.4.2 Spread option price 

The simplest spread option is an option of the European type on the difference of two 
underlying assets. The spread is naturally defined as the instrument S = {S(t)}t>o, whose 



n 



-rT 



ijj 



Cy^ j wkXk - k) 



k=l 
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value at time t is given by the difference 

S(t) = S 2 (t) - Si(t),t>0. (4.35) 
Buying such a spread is buying £2 and selling S\. 

The price of the simplest spread option is a particular case of (|4.34|) and equal to 

TV 

n = X \ X 2 ( 4 - 36 ) 

where 0jj denotes the European spread option price associated to the spread Y 2 — Y[. 

Vi,j = 1, . . . , N, Y{ and Y 2 J are log-normal underlying assets evolving according to the 
SDE (|4.15p . Let us denote the correlation coefficient between the two assets by p. It is 
possible to give a Black-Scholes type formula for the price of the European option with 
maturity T associated to the spread Y 2 — Y{ when the strike is K = 0, provided that the 
drifts (ii = (j>2 = r match the short interest rate r and the volatilities a\ and a 2 are constant 
in time. This formula was first derived by Margrabe in |41j as early as 1978, and it bears his 
name. It cannot be extended to the general case K ^ 0. 

Besides the fact that the case K = leads to a solution in closed form, it has also a practical 
appeal to the market participants. Indeed, it can be viewed as an option to exchange one 
asset for another at no additional cost. 

Proposition 5 When the strike K = 0, the European spread option price 0j j is given by 



4$(u4 j ) - ^>K) , (4-37) 



where 

M _ \n.{aP 2 /x\) 1 u f= jj _ ln(x J 2 /x\) 1 unf, 



and x\ = Yi(0), x ] 2 = Y 2 J (0), (<t* j ) 2 = (o-\) 2 — 2pa\a 2 + (a 2 ) 2 ,& is the standard normal 
cumulative distribution function, T the maturity, and uj = 1 for a call and u = —1 for a put. 

The proof is straightforward. 

4.4.3 Option on a weighted geometric average of a basket 

Let us consider an option of European type on the basket of securities of Eq. fj3.2j) with 
maturity T and strike K. The short-term interest rate is denoted by r and is assumed to be a 
deterministic constant. Then, if u = 1 for a call and cj = —1 for a put, the option value fj3.3j) 
can be written as 

n = e~ rT [ (u) [(x^ 1 • • • x^)»i+-"+»n - k] } + Ps(T)(x)dxi ■■■dx n (4.38) 
Jr™ L J J 
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where Ps(T) 15 the joint density of the random variables Si(T), ■ ■ ■ , S n (T) and is given by Eqs. 
HUSH - P~T9"]) . We have that 

N 

n= Yl X i 1 --- X n T k 1 ,-,k n (4.39) 

where T^ ... fc n denotes the European option price at initial time t = associated to the 
instrumental geometric-average basket (if 1 * • • • ! 

Particular case 

Let us now consider the particular case n = 2. The European option on weighted geometric 
average is then equal to 

N 

n = x \4 r*j (4.40) 

where I\ j denotes the European option price at initial time t = associated to the instru- 
mental geometric basket ( Y^ 1 Y 2 2 ' 



Recall that and Y^ , V«, j = 1, • • • , iV are lognormal underlying assets evolving according 
to the SDE (|4.15p . If the drifts fj,\ = \i 2 = r and the volatilities a\ and a 2 are constants 

in time, the price of the European Call option with maturity T associated to the basket 

i 

(YpYf 2 Y 1+W2 when if = is given by a closed form formula. 

Proposition 6 When the strike K = 0, the European Call option on a weighted geometric 
average price Tjj is given by 



Tij = exp(-rT)Y{(0) mWl Y 2 j (0) mW2 exp { Ur - \af)wi + (r - \a j 2 )w 2 

\ afwf + a 2 i2 w^ + 2pa\a 3 2 w 1 w 2 w 2 T^ 

where w = w + w and p denotes the correlation coefficient between Y± and Y 2 . 
Proof. To ease the proof notation we shall omit i, j indices. 

E | [Y^T^Y^T)™ 2 }^^ } 



zuT+ 



(4-41) 



-rT 

e 



(4.42) 



= e - rT Yi(0)™ wl Y 2 (0)™ W2 el (r ~^ Wl+ ( r ~^ W2 }™ T E [e 7Z ] 
where Z is a standard normal variable and 



7 = y [a 2 w\ + a\w\ + 2pa\a 2 w\w 2 \ w 2 T. 
The result follows by using E(e 7 ^) = e^ 2 / 2 . 
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Remark 7 The derivations ^4-^4^ an d j4-39 ) show that a dynamics leading to an n- dimensional 



density for the vector of underlying asset prices that is the convex combination of n- dimensional 
basic densities induces the same convex combination among the corresponding option prices. 
Furthermore, due to the linearity of the derivative operator, the same convex combination 
applies to option Greeks such as delta or gamma. 

Remark 8 The results of this section can be easily extended to hold in the case of shifted 
lognormal densities J2^. 

5 Markovian projection 

In this section, we introduce a model that we call "Multivariate Uncertain Volatility Model 
(MUVM)". We prove that the MVMD model is a multivariate markovian projection of the 
MUVM. 

5.1 Multivariate Uncertain Volatility Model 

This is a model specified through a system of SDE s of the form 

d£i(t) = m &(t)dt + aft) Zi(t)dWi(t), i = l,...,n, (5.1) 

where each Wi is a standard one dimensional Brownian motion, m are constants, a% are 
independent random variables representing uncertain volatilities. Each erj takes values in a 
set of iV deterministic functions a\ with probability A^ {a\ and A^ as defined in the previous 
section) and o~i(t) denotes its generic value. We thus have, for all times in (e, +oo), with small 
e, 

(t i — > o~l(t)) with probability A \ 
(t i — ► of (i)) with probability A^ 

(t i — ► CJi(t)) = < 

(t i — ► erf (t)) with probability Af 

We assume that all the above volatilities for asset i have a common time-path from to e/2, 
and then from the reached common value <Tj(e/2) at time e/2 each time- function connects to 
the relevant o~f(e) to continue then as a\. This is an initial regularization that is needed to 
make the dynamics smooth and ensure existence and uniqueness of solutions for the related 
equation. If e is small and the volatilities are smooth in time then we may neglect the initial 
regularization when computing expectations. We also assume that randomness of the time 
functions is realized at time e/2. Hence the uncertainty is quite short-lived and after that 
every asset follows a geometric Brownian motion. This makes MUVM a quite stylized model, 
and indeed MVMD, whose link with MUVM we are going to clarify now, is definitely more 
interesting. For an analogous analysis of the univariate case and a discussion see [5]. 
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We assume that the assets £j are pairwise correlated through the driving Brownian mo- 
tions. To be more specific we assume that d(Wi, Wj) = Pijdt. 

5.2 Markovian projection 



Lemma 9 (Gyongi's Lemma). Let us consider an n-dimensional stochastic process (£ 



starting from with the ltd form: 

^ t = ^it)dt + v{U t )dW t , (5.2) 

where W_is a standard d- dimensional Brownian motion, j3 is an n-dimensional bounded process 
and v is an n x d bounded process with vv T being uniform positive definite. There exists a 
Markovian n-dimensional process (X.t)t>0 which has the same distribution as (£ t )t>o at each 
fixed single time t, and which is a weak solution to the following stochastic differential equation: 

dX t = [±(t,X t )dt + a(t,X t )dW t ,X = 0, (5.3) 

where 

aa T (t,x) = E \vv T \f = x\ ,fi(t,x) = E \p\£ = x\ . (5.4) 



X is called the Markovian projection (in dimension n) of £. 

Theorem 10 The MVMD model is a markovian projection in dimension n of the MUVM. 

Proof. The system of SDE (|5.ip can be written in the following manner 

d£(t) = diagifj.) £(t) dt + diag(§(t)) A d\Z\, Z n ] T (5.5) 

with [Z\, Z n ] T a vector of n independent standard brownian motions and A the Cholesky 
decomposition of the correlation matrix X defined by Sjj := CTiCTj pij. 

The MVMD model given by flHD, flOQD and (I4T2T]) can be written as 

dS(t) =diag{y) S(t) dt + a{t,S(t)) d[Z u Z n ] T (5.6) 

where o~(t,S_(t)) := diag{S_{t)) C(t, 5(t))B, where B satisfies BB T = Rwith Rthe correlation 
matrix defined previously by Rj j := pij. 

The Markovian process S_ verifies: (i) S_ and ^ have identical one-dimensional distributions, 
i.e. they have identical distributions at every single time t conditional on the common initial 
condition at time 0. (ii) The following equality holds: 

E[vv T \£(t) =x] = a a T (t,x). (5.7) 
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Indeed, Denoting v(t,£(t)) = diag(£(t))A., we have 

E[dm<?(£(i)) S diag(£(t)) 
E[ot ICC*) = £l = ^77 i = ( 5 - 8 ) 

where 5 is a Dirac delta function. 

It is easy to see that the density of the multivariate process £ is 

N 

Pt(t)(2)= Yl A l 1 "- A n" P\ Y k Ht) ...Y kn (tw(£) ( 5 ' 9 ) 



'[Y 1 K1 (t),...,Y^(t)]T' 



and then we have after calculation 

N 



S p £(t) (x) = ^ A*..# V fcl '-^(t) P[ ^ 1(t)i>>>iyifn(t)]T (x) (5.10) 



^1 )**• j^Ti — 1 

where 

■yfel,...,fe„ 



(t) = ^(t)«i^(t) .. 1 (5.11) 

L J J i,j = l,...,n 



It follows from (JESJ) and (f57TO|) that 

AT 



E[%t)^}]= A^.-.A^" P [ yfi (t) ,„., Y) fn (t) ]T(2) ( 5 - 12 ) 



^1 3 • • -fen — 1 



and 



E[diag(0)) £ diag(£(t)) S{£(t)=x}] 



N 



fcl .... — 1 

(5.13) 

A substitution in (|5.8p gives 

E[ot t |^) = x] = a a T (t,x). 

We conclude by using Gyongi's Lemma [9] ■ 

Corollary 11 The process £ has the same distribution function as the Markovian process S_ 
for any time t. Then the MUVM can be used instead of the MVMD model to price European 
options. 

Remark 12 The MUVM features the following interesting properties: 

• explicit dynamics; 

• explicit density function; 
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• semi- analytic formulas for European-style derivatives. 

• semi- analytic formulas for early exercise derivatives (eg American Options). 

The last property follows via an iterated expectation, with the internal filtration referencing 
time e/2, and is not shared my the Markovian projection MVMD. See again for a discussion 
of the univariate case. 

6 Copula function in MVMD 

The aim of this section is to determine the copula function associated to the MVMD model 
given by Eqs. (|4.4|) . (|4.2Up and (|4.2ip where each evolve lognormally according to the 
SDE (l4TT5j) . 

To this end, let us begin by introducing Sklar's theorem in n-dimensions. This theorem 
elucidates the role that copulas play in the relationship between multivariate distribution 
functions and their univariate margins. 

Sklar's theorem 13 Let H be an n-dimensional distribution function with margins F\, F 2 , ■ ■ ■ , 
Then there exists an n-copula C such that for all x in R n , 

H(x 1 ,x 2 , ...,x n ) = C (Fi(x 1 ),F 2 (x 2 ), . . .,F n (x n )) . (6.1) 

If Fx, F 2 , . . . , F n are all continuous, then C is unique; otherwise, C is uniquely determined on 
RanF\ x RanF 2 x . . .RanF n , where RanF denotes the range of the function F. Conversely, 
if C is an n-copula and F±, F 2 , ■ ■ ■ , F n are distribution functions, then the function H defined 
by k6.1\) is an n-dimensional distribution function with margins F\, F 2 , . . . , F n . 

Proof. See Nelsen |43| for a detailed proof of this theorem. ■ 

Let us now define "quasi-inverses" of distribution functions 

Definition 14 Let F be a distribution function. Then a quasi-inverse of F is any function 
i^C -1 ) with domain [0, 1] such that 

1. if t is in RanF, then F^~^(t) is any number x inW such that F(x) = t, i.e, for all t in 
RanF, 

Wi^- 1 ^*)) =t 

2. if t is not in RanF, then 

= mf{x\F(x) >t} = sup{x\F(x) < t} . 

If F is strictly increasing, then it has but a single quasi-inverse, which is of course the ordinary 
inverse, for which we use the customary notation F . 
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The following corollary can be found in Nelsen |43| p. 47. 

Corollary 15 Let H, C,F 1 ,F 2 ,..., F n be as in Theorem\M, and let f[' X \ F 2 (_1) ,. . . , F^f 1} 
be quasi-inverses of F\, F 2 , . . . , F n , respectively. Then for any u in [0, l] n , 

C(u 1 ,u 2 ,...,u n ) = H^Ft 1) (n 1 ),F 2 i - 1 \u 2 ),...,Fj- 1 \u n )^ . (6.2) 



We are now able to determine the copula function associated to the MVMD model. 

Proposition 16 The copula function associated to the MVMD model, can be written as 

N 

C( Ul ,...,u n )= Yl A i 1 - A n n ^m(hi(F Si \ t) ( Ul )),...,h n (F^\ t) ( Un ))) (6.3) 

3>m denotes the standardized n- dimensional normal distribution function with correlation ma- 
trix M given by 

vHt)v^{t) 

for i,j G {l,...,n}, where V^{t) = <J jjj a^{s) 2 ds. is the inverse of the cumulative 

distribution function F s .^ of Si(t) given by : 

- r 1 / r i \i 

%W = E A * * [^y ( ln ^oj " ^ + 2^ (t)2 )J ' (6 ' 5) 

where $ is the usual standard normal cumulative distribution function, and hi is given by 

hAx) = (in— ^ ^t + -V^{t)A . (6.6) 

Proof. Using Corollary [15l we see that the MVMD copula is 

N 

C{ui,...,u n ) = Y A^...A^F [y fc 1(i)j ^ (6.7) 

h\ , . . . , kn — 1 

where F ^ t „ T is the cumulative distribution function of the vector (t), Y^ n (t)] T . 

Because i^t, ,,, ,,fc„,^ 1T is a cumulative distribution function and each Y- ki evolve log- 
[Y 1 1 (t),...,Y n n (t)] 1 i 

normally, it follows that 

V(xi,...,x n ) G R n , F [Y ^(t),...,Y^(t)]T (xi,...,x n ) = $ M (hi(xi), ...,h n (x n )) (6.8) 
from which the copula's expression follows. ■ 
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Corollary 17 The MVMD copula is a mixture of multivariate copulas that are the standard- 
ized multivariate normal distribution with correlation matrix M given by ( ffi.^| ) and margins 
G\ , . . . , G n defined as follows : 



Gi{x) = F m 



< x„ [ Vt{t) x + ln(lf *(0)) + - \vtit) 2 



(6.9) 



Proof. <]?m is the standardized n-dimensional normal distribution function with correlation 
matrix M given by (|6,4p . 

We prove now that each <I>m (h\(Fg^(ui)) , •••,h n (Fg 1 ^(u n )) \ is a copula. 

Because <I>m is an n- dimensional distribution function, we only need to prove that each 
hi o Fgf t \ is inverse of a univariate distribution function and then the result is deduced using 
Sklar's theorem. To this end fix an % € {1, • • • , n} and let Gi be the function from R to [0, 1] 
given by 



G i (x) = F Si{t) 



exp (v*(t) x + ln(^(0)) + ^t - \v^(t) 2 ^ 



Gi is a distribution function. Indeed 

1. Gi is increasing as composition of increasing functions; 

2. lim^-oo Gi(x) = Hnxr-K) F Si ( t ) (») = 0; 

3. lim^+oo Gi(x) = ]hn x ^ +00 F s .^(x) = 1. 
And its inverse G~ l = hi o F^ 1 ^. m 

Corollary 18 The MUVM has the same copula function as the MVMD model. 
Proof. This is a consequence of the markovian projection property. ■ 

7 Numerical Results for pricing: SCMD vs MVMD 

In this section we present some results for the pricing of three typical options : European 
Call on a weighted arithmetic average containing only long positions, European Call Spread 
option and European Call option on a weighted geometric average of a basket. We inves- 
tigate these options in the SCMD and MVMD frameworks in order to compare them. The 
performance of our approach is investigated by comparing the prices under the two models. 

For numerical sake, we focus on the two dimensional case where each individual component 
of the asset is modeled with a mixture of two lognormal densities. We assume also that the 
short-term interest rate r is deterministic and constant throughout the life of the option (i.e., 
until the maturity date T). Then, from Eq. f|3 . 3[) . the European Call prices tested in this 
section are given by the risk-neutral expectation 

n = e~ rT E [{B T - K) + ] (7.1) 

where B is the underlying basket instrument at maturity T. 
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7.1 European option on a weighted arithmetic average of a basket 

The European Call prices tested in this section are given by (|7.1|) where B is given by 
(|3.ip with (v)k)k=l,2 > f° r the option on a weighted arithmetic average containing only long 
positions. We call this option "Vanilla basket". Instead, B is given by Eq. (|4,35p for the 
spread option. 

Note that, under MVMD, the vanilla basket option price is given by Eq. f|4.34j) with 
n = N = 2 and the spread option price is given by Eq. (|4.36p with N = 2. 

The parameters of the test baskets are given in Table [TJ The interest rate r is 5%. The 
time to maturity (T) is one year. The strike K takes the values K = 0.7, K = 1 and K = 1.3. 
In order to obtain the fair price of the options under SCMD, 100,000 Monte Carlo runs are 
performed and an Euler scheme with time step At = 1/360 is applied. The first comparison 
uses a correlation p = 0.6. The results are given in Tabled The second comparison is done 
for a correlation p = 1. The results are shown in Tabled The standard error of the prices is 
given in the parenthesis. 





Vanilla Basket 


Spread 


Initial prices ([Si(0), 5 2 (0)]) 


[1,1] 


[0.7,1.7] 


drift ([m,^]) 


[5 %,5 %] 


[5 %,5 %] 




[0.6,0.4] 


[0.6,0.4] 


[AjiiA!] 


[0.7,0.3] 


[0.7,0.3] 


KV?] 


[0.3,0.2] 


[0.2,0.1] 




[0.25,0.35] 


[0.4,0.5] 


weights [w\,W2\ 


[0.5,0.5] 


[-1,1] 



Table 1: Basket option parameters 
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if = 0.7 




Vanilla Basket 


Spread 


MVMD 
SCMD 


0.3380 (0.0007) 
0.3386 (0.0007) 


0.4413 (0.0019) 
0.4365 (0.0019) 


K = 1 


MVMD 
SCMD 


0.1202 (0.0005) 
0.1200 (0.0005) 


0.2868 (0.0017) 
0.2833 (0.0017) 


K = 1.3 


MVMD 
SCMD 


0.0290 (0.0003) 
0.0296 (0.0003) 


0.1810 (0.0014) 
0.1836 (0.0014) 



Table 2: European Call on Basket Prices and Standard Errors for p = 0.6 



if = 0.7 




Vanilla Basket 


Spread 


MVMD 


0.3404 (0.0008) 


0.4199 (0.0018) 


SCMD 


0.3411 (0.0008) 


0.4193 (0.0019) 


K = 1 


MVMD 


0.1307 (0.0006) 


0.2611 (0.0016) 


SCMD 


0.1305 (0.0006) 


0.2647 (0.0016) 


K = 1.3 


MVMD 


0.0364 (0.0003) 


0.1661 (0.0013) 


SCMD 


0.0373 (0.0003) 


0.1637 (0.0013) 



Table 3: European Call on Basket Prices and Standard Errors for p = 1 

We shall prove, in a later paper [9], that Kendall's tau is different between the SCMD and 
MVMD models especially when p is high. Here we report the result without proof. 

Theorem 19 (Kendall's tau for MVMD). 



r (Si,S 2 ) 



2 

IT 



of arcsin(p) 



i=i 

4 4 



i=l 



+ 4 Yj Yj aia i Wp* ( mx v ' my ^ ) + ^ph ' ~ m ^j 

1=1 j=i+l 



where 



Vx % - (J>Xi 



my. 

1 XI 
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pij = 



with 



4 + 4 v 4 + 4, 



\ 1 \ 1 \ 1 \2 \2\1 \2\2 

"l = A 1 A 2 , «2 = A 1 A 2 , a 3 = A 1 A 2 , «4 = A X A 2 



Px x = Px 2 = ln(5i(0)) + (/a - ^(<rl) 2 J T , Px 3 = Px 4 = ln(Si(0)) + (/it - ^(4 2 

p Yl = PY 3 = ln(5 2 (0)) + L 2 - T, pr 2 = m = ln(5 2 (0)) + L 2 - \{o\f 

<*x x = cx 2 = o'i'yT, °Yi = <7y 3 = ^VT, ax a = &x 4 = ^fv^, 0Y 2 = <Ty 4 = cr|vT. 

We see, from the above theorem, that Kendall's tau is given by a closed form expression 
under the MVMD model. This is not the case for the SCMD model as we need simulations 
to get a numerical estimation of Kendall's tau from Monte Carlo. 

We show numerically in [9j that there is more terminal dependence in absolute value in 
SCMD than in MVMD. In the SCMD Kendall's tau does not change with the maturity, 
whereas, its absolute value goes down significantly as the maturity increases in the MVMD 
model. The relative difference of Kendall's tau between SCMD and MVMD is increasing with 
the maturity. It is relatively limited when p = ±0.6 and we could see more of difference when 
p=l. 

Despite this difference, the basket option price is not very sensitive to the difference 
between the two models, and indeed Table [3] shows that the prices obtained by the two models 
are close. Table [2] shows that this feature is maintained for a correlation p = 0.6. Notice that 
the prices obtained by the two models when dealing with a basket option with long positions 
are closer to each other than those obtained when dealing with a spread option. From all the 
experiments we made it seems that the price of the basket option with long positions increases 
as the correlation between the assets increases for all strikes whereas the price of the spread 
option decreases. So intuitively, a spread option will be more expensive if the correlation 
between the assets is lower, whereas a basket option with long positions will be cheaper. It 
can be seen that increasing the strike decreases dramatically the prices of the two options in 
the two models for the different values of correlation. The price of the spread option is higher 
than the price of the basket option with long positions and the difference between the two 
option prices becomes smaller as the correlation increases. These features hold for all strikes 
in the MVMD and SCMD models and are quite reasonable. 

These results seem to suggest that an option on a weighted arithmetic average of a basket 
containing only long positions and a spread option do not depend in an extreme way on 
dependence between the different assets since even models that give different Kendall's tau 
(See Brigo, Rapisarda and Sridi [9 J or Theorem [19] above) give quite similar prices. 
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Because the MVMD and SCMD models give similar numerical results in pricing European 
Call option on a weighted arithmetic average containing only long positions and European 
Call Spread option, the MVMD model is the most convenient here since it allows to compute 
the option price in one single Monte-Carlo step which can then be evaluated rapidly. 

In the next section, we will price an European Call option on a weighted geometric average 
under the SCMD and MVMD models and investigate if this option is more sensitive to the 
statistical dependence difference between the two models. 

7.2 European option on a weighted geometric average of a basket 

The European Call price tested in this paragraph is given by (|7.ip where B is given by 
Eq. ([3~2]) . Note that, under MVMD, this option price is given by Eq. ([ITlOj) with N = 2. 

Table U] reports the parameters we use. The interest rate r is 5%. The time to maturity 
(T) is one year. The strike K takes the values K = 0.7, K = 1 and K = 1.3. In order to 
obtain the fair price of the options 100,000 Monte Carlo runs are performed and an Euler 
scheme with time step At = 1/360 is applied. The first comparison uses a correlation p = 0.6. 
The results are given in Table [5j The second comparison is done for a correlation p = —0.6. 
The results are reported in Table [6j A last comparaison uses a correlation p = 1. The results 
are shown in Table [71 



Initial prices ([5i(0), 5 2 (0)]) 


[1,1] 


drift (\pi,p 2 \) 


[5 %,5 %} 


[Ai,Af] 


[0.6,0.4] 


[A2,Al] 


[0.7,0.3] 


KVfl 


[0.3,0.2] 




[0.25,0.35] 


weights [w\,W2\ 


[1,1] 



Table 4: Basket Option parameters 

While Kendall's tau is different between the SCMD and MVMD models especially when p 
is high (See Brigo, Rapisarda and Sridi |9] or Theorem 1 19p . the option price is not as sensitive. 
That is Table [7] shows that the prices obtained by the two models are close. Tables [5] and 
[6] show that this feature is maintained for a correlation p = ±0.6. We see that the prices 
obtained by the two models are close to those obtained previously in Section [7.11 when dealing 
with an option on a weighted arithmetic average of a basket with long positions. All the 
experiments show that the price of the option increases as the correlation between the assets 
increases for all strikes. It can be seen that increasing the strike decreases dramatically the 
prices in the two models for the different values of correlation. 
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K = 0.7 




Option price 


Standard error 


MVMD 


0.3313 


0.00074 


SCMD 


0.3312 


0.00075 


K = 1 


MVMD 


0.1154 


0.00055 


SCMD 


0.1159 


0.00057 


K = 1.3 


MVMD 


0.0267 


0.00028 


SCMD 


0.0268 


0.00029 



Table 5: European Call on Basket Prices and Standard Errors for p = 0.6 



if = 0.7 




Option price 


Standard error 


MVMD 


0.3049 


0.00037 


SCMD 


0.3045 


0.00037 


K = 1 


MVMD 


0.0584 


0.00025 


SCMD 


0.0574 


0.00025 


K = 1.3 


MVMD 


0.0016 


0.00003 


SCMD 


0.0013 


0.00003 



Table 6: European Call on Basket Prices and Standard Errors for p = —0.6 

These results seem to suggest that an option on a weighted geometric average of of a 
basket is not very sensitive to dependence between the different assets since even models that 
give different Kendall's tau (See Brigo, Rapisarda and Sridi |9j and Theorem [19] above) give 
quite similar prices. This is basically the same result we obtained for the arithmetic average 
basket in the previous section 

Because the above observations show that the MVMD and SCMD models give similar 
numerical results in European Call option on a weighted geometric average of a basket pricing, 
it is better to use the MVMD model allowing to compute the option price in one single Monte- 
Carlo step, which can be evaluated rapidly. 
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if = 0.7 




Option price 


Standard error 


MVMD 


0.3387 


0.00083 


SCMD 


0.3413 


0.00084 


K = 1 


MVMD 


0.1308 


0.00063 


SCMD 


0.1307 


0.00064 


K = 1.3 


MVMD 


0.0367 


0.00035 


SCMD 


0.0376 


0.00038 



Table 7: European Call on Basket Prices and Standard Errors for p = 1 

8 Conclusions and perspectives 

We illustrated how to extend in a conceptually simple fashion an asset price model, the 
so called (univariate and possibly shifted) lognormal mixture dynamics, that has been shown 
to reproduce well general implied volatility structures commonly observed on the market 
[6l [8j QUI l25| I42| . This model is formulated in the space of the so-called local volatility 
models. The extension aims at inferring an analytic expression for the local volatility of a 
multivariate security (such as e.g. a basket of underlying assets) that is consistent with (i) the 
individual dynamics of each component of the security as deduced by that security volatility 
smile and (ii) a given instantaneous correlation structure between different securities. 

A naive approach would consist in connecting univariate lognormal mixure dynamics mod- 
els for each asset through a instantaneous correlation connecting the Brownian motions driving 
different asset dynamics. We refer to this approach as simply correlated mixture dynamics, 
SCMD. 

However, we improve this approach by extending the mixture dynamics to the multivariate 
case in a more radical way, leading to the multi-variate mixture dynamics, MVMD, implying 
a multivariate mixture rather than single univariate mixtures patched together by Brownian 
correlations. While this is perfectly equivalent to SCMD for single assets, the main practical 
advantage of our MVMD extension is that our approach allows for a semi-analytic pricing 
of European style derivatives on the multivariate security in a way that takes into account 
the smile structures of the individual component securities and reduces computational time, 
while staying arbitrage free. 

This is a arbitrage-free dynamical model with a great potential for consistently modeling 
single assets and baskets (or indices) volatility smiles. 

We also showed that our approach performs remarkably well in terms of basket option 
pricing with a smile structure of implied volatilities. 
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We further introduced a multivariate uncertain volatility model, MUVM, which is not in 
the local volatility family but that provides the same european option prices and having the 
same copula function across assets as our MVMD model. 

It is well known that linear correlation is not enough to express the dependence between 
two random variables in an efficient way. That is why we are led to use measures such as 
Kendall's tau and Spearman's rho. While we introduce a theorem illustrating the difference 
in Kendall's tau for SCMD and MVMD (and hence MUVM). in a later paper [9] we establish 
closed form expressions for tau and rho under the MVMD approach. Such expressions do not 
hold for the traditional SCMD approach for which we need simulations to get a numerical 
estimation of such quantities. Unlike Kendall's tau, we find that Spearman's rho related 
to the MVMD approach is given by a large closed form formula that makes the numerical 
implementation difficult to handle. That is why, numerically, we focus only on Kendall's 
tau: we compute tau in the MVMD approach and compare it to its estimation in the SCMD 
approach with the same parameters. The results show that the dependence structure for the 
two approaches is different especially for high correlations. 

In this paper, we have tested our approach in baskets of 2 stocks with a smile structure 
of implied volatilities. Future extensions include (i) the testing of this approach in actual 
situations as swap rates derivatives within the LIBOR Market Model. Such an extension 
would allow computing in a quasi-analytical fashion the swap rates smile given the smiles in 
the individual caplets and an instantaneous correlation assumption, (ii) As for the pricing of 
European options, the hedging would be easy, accurate and rapid. (Hi) Our approach allows 
to compute European option prices in one single Monte-Carlo step which can then be quickly 
and accurately evaluated. However, for a path dependent payout on the asset, given that 
one needs to discretize the dynamics anyway, the traditional SCMD approach may be more 
convenient. 
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